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Based on the Liouville-von Neumann equation, we obtain a closed system of equations 
^ ' for the description of a qutrit or coupled qutrits in an arbitrary, time-dependent, external 

magnetic field. The dependence of the dynamics on the initial states and the magnetic 
field modulation is studied analytically and numerically. We compare the relative en- 
tanglement measure's dynamics in bi-qudits with permutation particle symmetry. We 
find the magnetic field modulation which retains the entanglement in the system of two 
coupled qutrits. Analytical formulae for the entanglement measures in finite chains from 
2 to 6 qutrits or 3 quartits are presented. 
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', 1. Introduction 

Multi-level quantum systems are studied extensively, since they have wide appli- 
^ cations. Some of the existing analytical results-'^^ for spin 1 are derived in terms of 

r> ■ a coherent vectoJ^. The class of exact solutions for a three-level system is given 

, in Ref. |3l The application of coupled multi-level systems in quantum devices is 

actively studied^. The study of these systems is topical in view of possible appli- 
cations for useful work in microscopic systemi^. Exact solutions for two uncoupled 
qutrits interacting with the vacuum are obtained in Ref. [6l For the case of qutrits 
interacting with a stochastic magnetic field, the exact solutions are obtained in 
Ref. [71 The exact solutions for coupled qudits in an alternating magnetic field, to 
our knowledge, have not yet been found. 

The entanglement in multi-particle coupled systems is an important resource for 
many problems in the quantum information science, but its quantitative value com- 
puting is difficult because of different types of entanglement. Multi-dimensional 
entangled states are interesting both for the study of the foundations of quantum 
mechanics and for the topicality of developing new protocols for quantum commu- 
nication. For example, it was shown that for maximally entangled states of two 
quantum systems, the qudits break the local realism stronger than the qubitfP, 
and the entangled qudits are less influenced by noise than the entangled qubits. 
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2 E. A. Ivanchenko 

Using entangled qutrits or qudits instead of qubits is more protective from inter- 
ception. From a practical point of view, it is clear that generating and saving the 
entanglement in a controlled manner is the primary problem for the realization of 
quantum computers. Maximally entangled states are best suited for the protocols 
of quantum teleportation and quantum cryptography. 

The entanglement and the symmetry are two basic notions of quantum mechanics. 
We study the dynamics of multipartite systems, which are invariant at any subsys- 
tem permutation. The aim of this work is finding exact solutions for the dynamics 
of coupled qudits interacting with an alternating magnetic field as well as the com- 
parative analysis of the entanglement measures in a finite chain of coupled qudits. 
The paper is organized as following. The Hamiltonian of the anisotropic qutrit in 
an arbitrary alternating magnetic field is described in Sec. II. Then the system of 
equations for the description of the qutrit dynamics is derived in the Bloch vector 
representation. We introduce a consistent magnetic field, which describes an entire 
class of field forms. In section III we find an analytical solution for the density 
matrix in the case of isotropic interaction. Analytical formulae, which describe the 
entanglement in finite spin chains of qutrits or quartits, are presented in Sec. IV. 
The results are demonstrated graphically in Sec. V at specific parameters. The 
conclusions are given in Sec. VI. 



2. Qutrit 

2.1. Qutrit Hamiltonian and Liouville-von Neumann equation 

We take the qutrit Hamiltonian (for the spin particle with s=l) in the space of 
one qutrit in the basis |1 >= (1,0,0), |0 >= (0,1,0), | - 1 >= (0,0,1), in an 
external magnetic field h — (hi, h2, h^) with anisotropy, in the form 

Hilt) = hiSi + h2S2 + hsSs + Q{Sl - ^^^i^S2s+ix2s+i) + diSf - Sj), (1) 

where hi, /12, ^3 are the Cartesian components of the external magnetic field in 
frequency unitsfwe assume h = 1, Bohr magneton /i^ = 1); Si, 6*2, S3 are the 
spin-1 matrices^ i?2s+ix2s+i is the unity matrix; Q, d are the anisotropy con- 
stants. When the constants Q, d are zeros, then the two Hamiltonian eigenvalues 
are symmetrically placed with respect to the zero l evel . 

There exist many useful baseW Allard and Rar^AR) formed the Hermitian 
basis Ca from the linear combinations of the irreducible tensor operators. This 
basis is normalized so that Si = Ci^x,S2 = Ci^y^S^ = Ci^z, irrespective of the 
spin quantum number s. It is convenient to construct the spin Hamiltonian for any 
spin. Hereinafter we use the Hermitian b asi^^. From the physical point of view, for 
important physical applications the basij™ is preferred. It is not necessary for 
the basis to be Hermitian since the results of the calculations are independent of 
the choice of base, but there is a significant advantage of the Hermitian basis. It is 
useful that the Liouville-von Neuman equation does not involve complex numbers 
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Entanglement dynamics in finite qudit chain 3 

and can be solved using real algebra. It makes numerical calculations faster and 
simplifies the interpretation of the equation system. The transition matrix deter- 
mines the coupling between the generalized Gell-Mann and (AH) Hermitian matrix 
bases. This coupling for qutrit is presented in Ref. 9. 

The qutrit dynamics in a magnetic field is described in the density matrix for- 
malism using the Liouville-von Neumann equation 

idtp^[H, pI p{t = Q) = pQ. (2) 

It is convenient to rewrite Eq. ([2]) presenting the density matrix p in the decomposi- 
tion with a full se of orthogonal Hermitian matrices Cq, (further the summation 
over the Greek indices will be from to 8 and over the Latin ones from 1 to 8) 

p — —^CaRa- (3) 
V6 

Since Tr Q = for 1 < i < 8, then from the condition Trp = Rq it follows 
that i?o = 1. And although the results are independent of the basis choice, in 
this basis the functions Ri = Tr p Ci have the concrete physical meanin The 
values Ri, i?2, ^3 are the polarization vector Cartesian components; R4 is the two- 
quantum coherence contribution in R2; R5 is the one-quantum anti-phase coherence 
contribution in R2; Rq is the contribution of the rotation between the phase and 
anti-phase one-quantum coherence; i?7 is the one-quantum anti-phase coherence 
contribution in Ri ; Rg is the two-quantum coherence contribution in Ri . 
The Liouville-von Neumann equation in terms of the functions Ri takes the form of 
a closed system of 8 real differential first-orde r equat ions. This system of equations 
can be written in a compact form as followin^lSIIll 

dtRi = eijihiRj, (4) 

where e^/ are the structure constants, hi — 2{hi, ^3, ^^^^^^ '-'1 ^"^^ Hamil- 
tonian components Eq. ([T]) in the basis Cq. 

2.2. The consistent field 

Let us consider the qutrit dynamics in an alternating field of the form 

h{t) = {ujicn{ujt\k), ujisii{ujt\k), uJodn{ujt\k)) , (5) 

where cn, sn, dn are the Jacobi elliptic functioniSl Such field modulation under the 
changing of the elliptic modulus k from to 1 describes the whole class of field forms 
from trigonometric^^ (cn(cji|0) — cosuit, sn(wt|0) = sinwt, dn(a;t|0) = 1 ) to the 
exponentially impulse ones (cn(wt|l) ~ sn(cjt|l) = thujt, dn{ujt\l) = 

The elliptic functions cn{ujt\k) and sn{ujt\k) have the real period while the 
function dn{ujt\k) has a period of half the duration. Here K is the full elliptic 
integral of the first kind"*^^. In other words, even though the field is periodic with a 
common real period but as we can see, the frequency of the longitudinal field 
amplitude modulation is twice as high as that of the transverse field. We call such 
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field consistent. 

Let us make use of tfie substitution p = a^^rai witli the diagonal matrix 
Qfi = diag (/, 1, Z""'^), where f{ujt\k) — cn(ajt|fc) + isn(ajt|fc). Then we obtain the 
equation for the matrix r in the form 

idtr = [aiHa^^ - iaidt{a^^),r]. (6) 

The equation for the matrix r without taking into account the anisotropy can be 
written as following 

idtr = [uJiSi + S dii{ujt\k)S3, r], r{t = 0) = po, S ~ ujq — uj. (7) 

-^t k = equ ation ([7]) describes the dynamics of the qutrit in a circularly polarized 
The exact solutions of this equation are known, and under certain 
initial conditions the explicit formulae are given in Ref. TS", At the exact resonance, 
w = Wo it is straightforward to present ^ in the deformed field {k 0) ^ for the 
given initial condition p — po: 

p[t) = a-^ e Voe 'ai- (8) 

Explicit solutions for specific initial conditions are given in Ref. |9l 

3. Bi-qutrit 

In the space (E) the bi-qutrit density matrix can be written in the Bloch 
representation 

Q = ^Ra/sCa <X) CfS, RoO = 1, Q{t = 0) = Qq, (9) 

where (8) denotes the direct product. The functions RmO,Rom characterize the in- 
dividual qutrits and functions Rmn characterize their correlations. 
Let us consider the Hamiltonian of the system of two qutrits with anisotropic and 
exchange interaction in a magnetic field in the following form 
. ^. . ^ 

H2=H{h)(E) E2S+1X2S+1 + E2S+1X2S+1 (E) H { h ) + J S^ (g) S",, (10) 

where h and h are the magnetic field vectors in frequency units, which op- 
erate on the first and the second qutrits respectively, and J is the constant 
of isotropic exchange interaction. We study the dynamics of two qutrits in the 
consistent magnetic field h = {ijjicn.{uit\k)), a;isn(a;t|/c), uJodii{ojt\k)) , h = 
{wicn(v[7t\k), zuisn{u}t\k) , zuodn{ujt\k)) at the anisotropy constants equal to 0. 
Let us transform the matrix density g = a^^r2a2 with the matrix a2 — oti (E ai. 
The equation for the matrix r2 takes the form idtr2 — [i?(dn(tjt|fc)), with the 
transformed Hamiltonian H {dii{ujt\k)^. 

Since d'n{ujt\k)\k=o — 1, then the transformed Hamiltonian H docs not depend on 
time, and the solution for the density matrix in the circularly polarized field has 
the form 

0{t)^a^^e-'^'Qoe'"'a2\k=o. (11) 
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In the consistent field at resonance u — vjq = ujo = h at equal zui — uji the 
Hamiltonian eigenvalues equal to — 2J. —J, J, J— 2a;i, — J — wi, J — cji, —J + uji, J + 
wi , J + 2u;i . This allows to find the exact solution in the closed form for any initial 
condition since the matrix exponent e*^' in this case can be calculated analytically. 
For a larger number of the qudits with a pairwise isotropic interaction, the general- 
ization is evident. In the case of interaction of qudits with a different dimensionality, 
the reduction of the original system to the system with constant coefficients can be 
done by choosing, for example, the transformation matrix for spin-3/2 and spin-2 
in the form 

diag(/3/^ /-3/2)^diag(/^ /, 1, f-\ r'y m 

However, the Hamiltonian eigenvalues cannot be found in a simple analytical form 
because of the lowering of the system's symmetry. 

4. Analytical formulae for entanglement measures 
4.1. Entanglement in the bi-qutrit 

For the initial maximally entangled state which is symmetrical at the particle per- 
mutation 

IV' >= -/^ E I* > >' (13) 

^ z— — 1 

in the consistent field at the resonance uj = vjq = luq = h at equal wi — wi, 
the exact solution for the correlation functions is given in Ref. |21 The correlation 
functions have the property i?^^ = Rpa^ i.e. the symmetry is conserved during the 
evolution, since the initial state and the Hamiltonian are symmetric with respect 
to the particle permutation. 

Given the exact solution, one can find the negative eigenvalues of the partly trans- 
posed matrix = ® E^q (here T denotes the transposition): = — 
— ^VBQ + 28 cos 3 Jt — 16 cos 6 Ji, £3 = — ^ (5 -I- 4 cos 3 Jt) . The absolute value of 
the sum of these eigenvalues 



mvw 

i=l 

defines the entanglement measure (negativity) between the qutrit^^QI 

The en tang lement between the qudits can be described quantitatively with the 

measure^ 



msM = JJ~J^^i3 ^ RioRo]Y, (15) 

where D is the basis dimension (for qutrit D = %). This measure equals to for the 
separable state and to 1 for the maximally entangled state, and it is applicable for 
both pure and mixed states. 
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6 E. A. Ivanchenko 

That is why for the maximahy entangled initial state of two qutrits, the entangle- 
ment in the consistent field is defined by the formulae with the found solution for 
the density matrix 

msM = ^\/4457 + 2776 cos 3 Jt - 632 cos 6Jt - 56 cos 9Jt + 16 cos 12 Jt (16) 
81 

This measure is numerically equivalent to the measure TO\/vi/l2ni22J -^j^jc}! defined 
by the absolute value of the sum of the negative eigenvalues. 

According to the definition for iV-qudit pure state^, the entanglement measure 
equals to 



N 

i 

Vn 



N 

where Si — ~Tt pt logf, pi is the reduced von Neumann entropy, the index i numer- 
ates the particles, i.e. the other particles are traced out. We use the logarithm to 
the base b to ensure that the maximal measure is normalized to 1. The base b equals 
to 3 in the qutrit case. 

Since the qutrit reduced matrix eigenvalues equal to Ai = A2 = 1^^(5-1- 
4cos3Jt), A3 = ^(17 — 8cos3Jt), then the entanglement measure in the bi-qutrit 
takes the form 

3 

r,2 =-^Aaog3A,. (18) 

Normalized to unity the measure I-concurrence which is easy to calculate is defined 
by the formulaJ 



^i = \j2 — ^\/(l-Trp?) = i\/57-|-32cos3Jf-8cos6Ji, (19) 

where d = 3 for a qutrit, pi — -^CaRaO is the reduced qutrit matrix. 
The measures myw, rasu, '72, rai do not depend on the parameters of the con- 
sistent field, the sign of the exchange constant at zero anisotropy parameters. It 
should be noted that the Wootters entanglement measure (the concurrence) in the 
system of two qubits with an isotropic interaction in a circularly polarized field at 
resonance is also independent of the alternating field amplitude^^, but depends on 
the exchange constant J and the initial conditions only. 

The numerical solution of the Liouville-von Neumann equation shows that if an 
identical external field operates on every qudit, the free Hamiltonian and the inter- 
action Hamiltonian are commutative operators, the measures considered are deter- 
mined only by the symmetric two-body interaction with the interaction constant 
of J. If a different field operates on every qudit, there arises a broken permutation 
symmetry of the total Hamiltonian, which changes the entanglement dynamics. 
Thus it is possible to control entanglement by changing the parameters of an ex- 
ternal field. 
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Entanglement dynamics in finite qudit chain 7 

It is possible to show that the distance measur 

— PoY depends 

on the parameters of the consistent field. 

At a zero external field the entanglement measure (fT5|) takes the analytical form 
at equal non-zero anisotropy parameters Q = d = d = Q 



4 

^gfecos (fcv/9J2T8QjTT6Q2t), (20) 

fe=0 

where q„ = 44577^ + 11616QJ^ + 47392gV + SSSSSQV^ + 163072(34 j4 + 
194560g5 J3 + 221184Q6J2 + 131072QV + 65536Q^ 

= 8 J2( J + 2Q)2 (347 J4 + 518Q J3 + 1440QV2 + 1504Q'V + 1024Q4) ; 
q2 = -8J2(J + 2Q)2 (79 J4 + 76Q J3 + 320Q2 + 448Q3J + 256Q''); 
93 = -8J3(7J - 4Q)( J + 2Q)3( J + 4Q), q^ = 16 J4( J + 2Q)4. 

4.2. Entanglement in the chain of qutrits 

We consider the Hamiltonian of the chain of N qutrits with the pairwise isotropic 
interaction in the consistent field h{t) ([5]) at resonance in the following form 

N-l N-2 

Hn = Y{h{t)t ®E®---®E + ®E®---®E), (21) 

where the summation is over different possible positions of in the direct products. 
Because the maximally entangled state of N qutrits 

1^ >^= 7f E N (22) 

and the Hamiltonian (PT|) have a permutation symmetry, it follows that the density 
matrix of N qutrits has symmetric correlation functions. 

The entanglement measures for many-particle multi-level quantum systems have 
not been studied enough and are difficult to calculate in the analytical form, that 
is why we will present analytical formulae only for the entropy measure t^n, which 
is defined by the eigenvalues of the reduced one-particle matrices for each qutrit. 
As the result of the mentioned symmetry the reduced matrices are equal to each 
other. Therefore the entanglement measure for N qutrits reads 

3 

r^AT = -^r^loggri. (23) 

i=l 

The eigenvalues of the reduced matrices for 3, 4, 5, and 6 qutrits are presented in 



{9J2 + SQJ+16Q^f\ 
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the table below 

N \ Ti ri = r2 7-3 



29-4 cos 5 Jt 17+8 cos 5 Jt 

75 75 

905-98 cos 3Jt-72 cos 7Jt 395+196 cos 3 Jt+144 cos 7 Jt 

2205 2205 

16919-1944 cos 5 Jt-800 cos 9 Jt 8687+3888 cos 5 Jt+1600 cos 9 Jt 

42525 42525 
21977-1694 cos 3 Jt- 1936 cos 7Jt-560 cos 11 Jt 9407+3388 cos 3Jt+3872 cos 7Jt+1120 cos 11 Jt 

53361 53361 



(24) 

The measures 773, 774, 775, 775 do not depend on the sign of the exchange constant 
like the measure 772- 

4.3. Entanglement in the bi-quartit 

The applied approach for qutrits is translated to qudits. For a spin-3/2 particle 
or a four-level system, also denoted as a quartit, we take the Hamiltonian in 
the space in the basis |3/2 >= (1,0,0,0), |l/2 >= (0,1,0,0), | - 1/2 >= 
(0, 0, 1, 0), I — 3/2 >= (0, 0, 0, 1) and use the matrix representation of a complete 
set of the Hermitian orthogonal operators. 

We will find the analytical formulae in the bi-quartit and in the 3 quartits (in 
bi-pentit, see below) with a pairwise isotropic interaction of the initial maximally 
entangled state in a consistent magnetic field at resonance without taking into 
account the anisotropy. 

The negative eigenvalues of the partly transposed matrix g^** are 
equal to Ai — TM^"-*^"^ ~ 12cos5Jt), A2 = A3 = A4 = A5 = 
- jig V409 + 288 cos bJt-72 cos 10 Jt, Ae = tm(~37+ 12cos5Ji). The entangle- 
ment measure in the bi-quartit equals 

^yiy"^*=ElA^I- (25) 

1=1 

The entanglement between the quartits is described quantitatively with the 

'21 
measure 

b^-qrt _ V1803365 + 191616 cos5 - 35808 cos 10 Jt - 6912 cos 15 Jt + 864 cos 20 Jt 
"^^^^ ^ 625V5 

(26) 

Since the quartit reduced matrix eigenvalues equal to Ai = A2 = ^^(13 -|- 
12 cos 5 Jt), A3 = A4 = -^{'^'^ ^ 12 cos 5 Jt), hence the measure 772 reads 

7;^«'-* = -VA,log4A,. (27) 



=1 



The I-concurrence is equal to 

bi—qrt 



1 



my ' = — %/553 H- 96 cos 5 Jt - 24 cos 10 Jt. (28) 
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Entanglement dynamics in finite qudit chain 9 

4.4. Entanglement in the quartit chain 

The eigenvalues of the reduced matrices for 3 quartits are equal to ri ^ r2 = 
0.141 + 0.068 cos ^ +0.04 cos 8 Ji, rg = r4 = 0.359 - 0.068 cos ^ - 0.04 cos 8 Jt. 
Therefore the entanglement measure for 3 quartits is following 

4 

= -Xl^'^°S4''i- (29) 



4.5. Entanglement in the bi-pentit 

For the spin-2 particle or a 5-level system, also denoted as a pentit, we take 
the Hamiltonian in the space in the basis |2 (1,0,0,0,0), |1 > = 

(0,1,0,0,0), |0>= (0,0,1,0,0), I -1>= (0,0,0,1,0), I -2 (0,0,0,0,1). 

The entanglement between the pentits is described using the measure^ 
mg^P"* = (0.802 + 0.106cos3Ji- 0.019 cos 4 Jt + 0.242 cos 7 Ji- 0.098 cos 10 - 
0.088 cos 14Jt + 0.067cosl7Jt-0.014cos 20 Ji)i/2. The I-concurrence is determined 
by the formulae = (0.791 + 0.114cos3Jt - 0.018cos4Ji - 0.005cos6Ji + 

0.230 cos 7Jt - 0.079 cos 10 Jt - 0.079 cos UJt + 0.060 cos 17 Jt - 0.015 cos 20 Jt)^/'^ . 
We have replaced the exact bulky rational coefficients by its decimal approxima- 
tions and the terms less than 0.001 have removed for inconvenience reduction. 
The pentit reduced matrix eigenvalues are equal to pi = P2 = gT25^^^'^'^ ~ 
140cos3Jt + 640cos7Jt - 448coslOJt), P3 = _P4 = gj^(513 + 280cos3Ji + 
320cos7Ji + 112coslOJt), P5 = -^(2753 - 280cos3Jt - 1920cos7Ji + 
672 cos 10 Jt), hence the measure r/2 reads 

^-EftlogsP. (30) 

i=l 

All the measures do not depend on the sign of the exchange constant and the 
parameters of the consistent field at zero anisotropy parameters. 



5. Numerical results 

Although the analytical expressions for the measures in a bi-qutrit myw, "niSM are 

different, but the numerical values are practically identical. The maximal deviation 

in the rectangle (1 > J > 0.01) x (100 >t>Q) equals 0.014. 

Measures rj2 and m/ qualitatively coincide with the measures mvw, "nisAi- 

We have found that the anisotropy of the qutrits disentangles them, namely the 

entanglement is decreased down to O.OOl^see graphs 1 and 2 in FigH]). 

In the constant longitudinal field h — — h = (0, 0, uiq) (the bi-qutrit Hamiltonian 
eigenvalues are equal to J, J, xi,X2, x^, — p, —p,p,p, where xi,X2,x^ are the roots of 
the equation x^ + 2x^J — p^x — 2J^ = 0, p = \/ ,P + Wg) Hamiltonian contains 
the asymmetric part, thus it follows that the density matrix for the initial symmetric 
state will not be symmetric because of the breaking of the symmetry of the particle 
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Fig. 1. Dynamics in the bi-qutrit: in the Fig. 2. Disentanglement of the maximally 

zero external field at equal anisotropy con- entangled state in the chain of 2,3,4,5,6 

stants Q = 0.0250, J = —0.1 (curve 1) and qutrits with J = 0.1. 

J = 0.1 (curve 2); in the consistent field at 

J = 0.1, Q = the curve 3 shows complete 

coincidence of myw ^'Hd msM', the curve 4 

is the measure 772; "i/ is the curve 5. 



permutations. The analytical solution is cumbersome. In the constant longitudinal 
impulse field it = -h = {0, 0, 2 (6'((t - 17){t - 60)) + 6'((40 - ^(57 '~t){t~ 60)))) 
the entanglement dynamics is blocked^ at wo 3> J . This points to the possibility 
to control the entanglement. 

In Fig|2] we present the comparative dynamics of the entropy measure in the finite 
qutrit chain. The disentanglement dynamics of the measures t?4, ?75, % is similar 
to the one in the case of two qutrits, but with smaller oscillation amplitude, i.e. 
larger number of the qutrits disentangles less than two qutrits (0.889 < 773 < 1). 

The measures in bi-quartit, as shown in Fig [31 qualitative coincide, almost com- 
pletely Tj2~'^^^ and my'^^* . The disentanglement 3 quartits is insignificant less than 
in bi-quartit. 

The disentanglement measures in bi-pentit, as shown in Fig|4l qualitative coin- 

•j 1 J. 1 J- 1 bi — pnt 1 bi—pnt 

cide, almost completely rrij and nigj^f . 



6. Conclusion 

The comparative analysis of the bi-qutrit entanglement measures on the base of the 
analytical solution for the density matrix demonstrates that, in spite of the different 
approaches to the derivation of the formulae for the entanglement, all the formulae 
yield quite close results (Fig. [T|), and the measures myw and msM are practically 
equal. This is in accordance with the general results for the entanglement in the 
systems with a permutational symmetrji^. 

The analytical formulae for the measures ^73, 774, 7?5, rjQ are similar to the measure for 
two qutrits 772, but with a numerically smaller oscillation amplitude, i.e. the larger 
number of the qutrits disentangles fewer than two qutrits. 

Nevertheless, the comparison of measures in two coupled qutrits, quartits, and 
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Fig. 3. Disentanglement of the maximally 
entangled state in the chain of 2, 3 quar- 
tits with J = 0.1. The measures m'j~'^^* , 

"*Sm''''*' ^2*"''^*' "*v'vi^''*' presented by 
the curves 1, 2, 3, 4 respectively; r)^* is the 
curve 5. 



Fig. 4. Disentanglement of the maximally 
entangled state in the chain of 2 pentits 
with J = 0.1. The measures mj' 

™Sm''"*' '?2'~'"''' ^""^ presented by the 
curves 1, 2, 3 respectively. 



pentits 

hi — qutrit mvw — rnsM 

bi - quartit r?^'""'* ^ 771^^^""' (31) 
bi - pentit m%^* ^ mj*"^"* 

on the base of analytical solutions shows the absence of a full coincidence of the 

measures even in a particular case of disentangling a maximally entangled state. 
In other words, it is impossible to prefer any measure, there remains therefore the 
question concerning the quantitative determination of entanglement even in case of 
two multi-level particles. 

The author is grateful to A. A. Zippa constant invaluable support. 
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